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Program Title: Multiplication of Matrices (Relocatable) 

1. PURPOSE 

1.1 To compute in Floating point form the product matrix 

C = AB 

1.2 To compute in floating point form the product matrix 

C = A*B 



2. RESTRICTIONS 

2.1 Maximum matrix size approx, 17 x 17. May be changed by varying address 
constants (see 4,5) 

2.2 Routine dependent on matrices being stored in memory prior to operation 

Matrix A - stored rov;-wise 02 00- on 
Matrix B - stored row -wise 1400- on 

2.3 No output - Product matrix is stored row-wise 4000-dn 

3 . METHOD 

3.1 AB=:C C. . = 1^^ A,^^ B^^. 

^(mxn) ~ (mxl) (Ixn) 

3.2 A*B=C C. . = ^ . 

k-1 kx kj 

C = (A ) B 

(l2ai) ~ (Ixra) (mxn) 



4. USAGE 

4,1 m, n, 1 must be entered in numeric format in locations 0^35, 0136, 0137 - 
binary 39 



4,2 Options 

The proj 

start at 0013,0, and it will automatically make the necessary changes. 



4.2,1 The program is set for finding AB. If you want to change it to find A B, 



4.2.2 To change the program to find C = AB once it has been set to find C = A B, 
start at 0033,0, 

4.2.3 If you don't need to change the program, start at 0037,0. 
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4.2.4 The routine stops on an HTR 0037.0 in location 0113.0. 

This nay be replaced by any other exit if desired. 

4.3 Sl^orage. 

( 0000 - 0017 - Routine to change program from C=AB to C=A B 

< t 

( 0020 - 0036 - Routine to change program from C=A B to C=AB 
Relocatable ( 

( 0037 - 0121 - Program 

( 

( 0130 - 0161 - Constants and modified addresses 

Remaining storage free for matrices or other uses, (see 2,2, 2.3, 4.5). 
Note: Relocation does not affect matrix addresses. 

4.4 Timing examples (either option) 

4. 4.1 (16 X 10) (10 X 16) =6-6x1^ 20 minutes 

4.4.2 (10 X 16) (16 X 10) = ^0 x 10) 15 minutes 

4.5 Routine may be modified to act on matrices in different locations by 
changing the following address constants: 

01.54 - (Matrix A address) 1/2 + 00 0000,0 00 xxxx.O 

01.55 - Matrix B address + 00 0000.0 00 xxxx.O 

01 56 - (Matrix C address) 1/2 ■*- 00 0000.0 00 xxxx.O 

For example, if you want A. A = c, you would store: 

01 54 * 00 00000 00 0100.0 ) 

) Matrix A address = 0200.0 
01 55 + 00 0000.0 00 0200.0 ) 

01 56 + 00 0000.0 00 2000.0 Matrix C address = 4000.0 

4.6 Both L and V loops are used by the routine. 

4.7 To restart for new multiplication, reset new m, n, 1, check your 
address constants (if you changed them) and start as in no, 4.2, 



000070 ~4-6"6V005D"."0 - ~ UTV ■ 
4-00,0n2.0 CLA 
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0001 .0 +11 .01 4l .0 MPY 
" +113.6660.6 XAR 



660 2.0 "+6 3". 1 MV6 ~ " SUB 
+^.01 61.0 STO 



_Q00J.Q_+11,gi3T.P___MPY 
+^3. 0000.0 XAR 

"000^ .6 " +66 ."6TV5 ". 6" ~ "stO" 

+OO.On^.O CLA 



QQQ5,0 +1 1.01^1 .0 MPY 
W3.6606.6 XAR 



oooS. +63 .61 X f ". 6" " 'sua' 
+m:oi4^.o ADD 



0007^.0 . +p L.aL5,^^P_. „ ,App. . 
+41.0025.0 ALS 



001 0.0 +03.orj'p.o sub 

+52.7752.1 TRP 



0011.0 +60.0133.0 STO 
+03 ."Ol 37.-O" " -WB" 



0012.0 +UU.UUU0.0 — =^=^ 
-00.0000.0 — 



0013.0 +6H.0000.0 _ CTL 

+ff5:oo^o:o "CFL" 



0014.0 +00.0010.0 clX 
..._+^&..0Q6j_.q_._stp_ 



.00 1 5.0 ,.+QQ,cQQ.lJ_*.0.„ _ _CL A_ 

+66.0101.0 STO 



0016.0 +57.0037.0 TRA" 
+00.0060.0 CLA 



0017.0 +00. 0000. 0_ — " 
-60.6666.6 — 



662670" ~+56^.~665076 " " cTv 

+Qp_._0 1 33_._0 _ __ CLA 



0021 .0 +11.01^1 .0 MPY 
+43.0000.0 XAR" 

0022. 6 +63.61 41 .0 ' ~ SIJ3 
+IT.OI37.O MPY 



0023.0 +43.0000.0 XAR 
^ +66.0145.0 STO 



+ 11.01_41_*q MPY ^_^R_^_ 

0025.0 +H^. 000 0.0 XAR 

+03.01^1.0 SUB 



0026.0 +60.0161 .0 STO 
.+04-^0 U-^^- - ADD. 

0027.0 4- 01 , 015 ^. — ^DD 



+1^1.0025.0 ALS 



0030.0 +03.0137.0 SUB 
+^2.11bP..] TRP 



.0 +60.0133.0 STO 
" +05. 01 35.0 SU3" 



0032VO +oor;ooD'o.o- ■---■ 

-00.0000.0 — 



.0 +6^.0020.0 CTL 
" +65.00'HU.O' -CFL" 

003ii-'. ~+o"a;oo3u:o- re A- 

+60. 0061.0 STO 



0035.0 +00.0031 .0 CLA 
" +^7". 001 5 ." T TR A- 



"OO36 . " +00 .0000 ; O - - - ■ 

-00.0000.0 — 



0037.0 +6^.00to.O CTL 

+57.776o.O'TRA- 



"ooTo.o +66.0050.0 CTV — 

±00.0133.0 _CLA_ 

OO^H .0 +11 .0\h] .0 MPY 
+I|3.00C)6.6 XAR 



0052.b +03.0141 ,0 SU3" 

+iuoj37,o MPr 

00^3.0 +43.0000.0 XAR 

+60.01 U5.O STO" 



.0 +00.0132.0 CLA 
+ 1 1.01^1.0 MPY 



00^5*0 +43.0000.0 XAR 
+03.0l4l .0 " SUB 



7o +6o.or6i"".o sTO" 

+01 .0145.0 . ADD 

0047.0 +01 .0154.0 ADD 
+4l .OO25VO TOTS 



I . +P0JD134^ CLA 

05K0 +n >OlHl >0 MPY 

+h^.0000.0 XAR 



005 2 . +03.0 f M 76 ' ' SllB " 
, . _+g)..0I5J_._0. _ STO. 



0053. +0 0. Q]6l >Q CL A 
+11 .0136,0 MPY 



005^.0 +H3. 0000.0 XAR 

+D1.0J_51..0-__AD_D_ 



0055.0 +57,01 PP. TRA 

+m57oi467i STA 



0056.0 +57,01^^6,0 TRA 

+jQQ.0J32,0-._CLA_ 



0057.0 +6H,0060.Q CTL 



+57,7760,0 TRA 



6o5oVo"+QT.OTto:o ~"add~ 

+60.0132,0 STO 



0061 ,0 +03,0137.0 SUB 
+52-; 7762.- r - - TRP- 

0062 .-0 ~+5T,"oa37,-0"- - tra- 

+50.0037.0 TRZ 



0063.0 +00,01^0,0 CLA 
+60.0132,"0~""ST0- 



:)o6i^.o +^6.oo7"'a,0'' XTv" 

+ 00.0133.0 CLA 



00(95.0 +1 1 ^01 4l .0 MPY 
+^•3. 6000,0 XAR 



)o66.o +6~3"."6i^i .'0 ~'SU3' 
+ii.on6.o MPY 



0067.0 +i^3.0000,0_ XAR 

+66.6r52,6 STO 



0070,0 +00.013^.0 CLA 
. __:.._.+! J_JD141_,_0__ MPY. 

0071 - +k3. OOPO ._0 _ _ XAR . 

+03.oi4i,o SUB 



0072.0 +01 ,0152,0 ADD 
+01 ,01 56,0. _ ADD 

0073.0 +i|l .0025,0 ALS 
" "+fe,0r53V0'"5TA 



007^.0 +30.0142.0 FCA 

+5IaOI51JD__tra PAGE_6. 



0075.0 +00.01^0.0 CLA 
" +60VCr32V0""BT0- 



0076.0 +30.01 30.0 rCTT" 

_+35^Qlte^Q_.fST_. 



.0077-ja-+6iUDlOCUO--ilIL. 
+57.7760.0 TRA 



'oibb76~+o6~6T33^o"''cLA' 

_ tPl.pj 40.0 _ ADD 

0101 .0 +60.01 33>0 STO 
+03.0155.0 SUB 

f 02" ;o "+52-.7763:o ■ - trf 

t5_It.OO3_I.O_ TRA 



0103.0 +50.0037.0 TRZ 
+60.01 W.O CLA 



f oi^;i5 "+So ;oT 33 .'0' ' 'STO 

_±00_._0J3J(^.0___CLA_ 

01 05.0 +Q1 .Ol^p.O ADD 
+6U.UI S4.0 — STO- 
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"01 06.0 +6"6;oi Tu;o~ ' rrv" 
__ .j:03.oj36.o__^UB 

0107.0 +52.0110.0 TRP 
+57.003 /.U TRTT 



01 TO :0" +50:0t3'37~. 0" ""TRZ 
+57.0111.0 TRA 



0111 .0 +30.0130.0 FCA 
+35.01'42.0 F'St 

oil 2 ;q--:^qo ;ol Ifo :6~~ ~CL A" 

+6o.Qn4.o STO 



011.3,0 +77*0037.0 HTR 
+00.0000.0 CLA 

'ori¥.o' +oo^oooo;o~ "---" 

-00.0000.0 — 



0115.0 +00.0000.0 

-00 ;oooo ;o" " "-'-- 

'] 1 6.0 +00 . oooo";o " ~ - - - 
-00.0000.0 — 



0117.0 +00.0000.0 — 

" -00.0000.0' ^^-^ 



0120.0 +11 .01 4l .0 MPY 

_-t!}3j:}0QO^„.XAR . PAGE 7 

0121 .0+01 .0155.0 ADD 

+57-77T5.1 ^^' ' '" " ' 

"01 22.0 +00.0000.0 '^^^^^^ 

-00 .0000 .0 - -- _ _ _ 

0123.0 +00 . 0000. ., . - - - , , 

-00.0600.0 -— 



012^1.0 +00.0000.0 - — 

..00.00aQ*0...-„-rr.. 

.01 25.O-+QQ*0D00*0_...rr.-_ 
-00.0000.0 



0126.0 +00.0000.0 

-QO .00.00. 0..,--:.-. 

.01 27.0 +0Q..OP.OO.O .— -. 

-00.0000.0 — 



0130.0 +00. 0000; 
-oo.goop.p_ 

0131 .0 +00.0000.0 
-UU. 0000.0 



0132.0 +00.0000.0 
-00.0000.1 

0133.0 +00.0000.0 
— — ~ "-TOO. 0000.1 

013^.0+00.0000.0 

-00.0000.1 
013^.0 +00.0000.0 



-00.0000.0 

01 36. 0+00. 0006. 6 

-QQ.00Q0.0.. 

0137.0 +00.0000.0 — 



-00.0000.0 — 



1 \Q) .0 +00 ."0000 . 
=:a:).Qoooj 



Oli^l .0 +00.0000.0 
-00.0001 .0 

01^^2.0 +00.0000.0 
-OQ ,OQQQ.O 



01^3.0 +00,0000.0 
-00.0000.0 

oi'if^ .0 +oo.oooD-;o" 

-00.0000.0 



ill45^_Q +Qp_._00gp.0_ _ _ _ |*AG^_5 

"-oo.oooo.'o — ■ 
01I1S.0 +30.0200.0 rcA ' '" 

+07.1^00.0 FMP 

01^7,0 +0!l. 01 ^2.0 FAD 

+35.61^2.6 FST 



0150.0 +^7. Ill 6.] TRA 

015UO +00.0000.0 — 
-00. 0000. ~— 



0152.0 +00.0000.0 ~— 
-00 ,0000 . . - -_-_ 

0153.0 +•^.1+000.0 FST 



0154.0 +00.0000.0 ^^^ 
-00.0100.0 — 



0155.0 +00.0000.0 — 



0156.0 +oo.o6uu.u =^^ 
0157.0 +00.0000.0 — 

"" ~~"-ua:oouo:o"-"---- 



01 6o;o"+o'o:oooo-.cr~----~ 
-00.0000.0 .— 



01 61 .0 +00.0000.0 — 
-00.0000.0 ::^^^ 



